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Bayes  and  Equivariant  Estimators  of  the  Variance  of  a Finite  Ibpulation 


By 

S.  Zacks  and  H.  Solomon 


0.  Introduction. 

Let  x.^, ...  ,x^  be  the  values  of  a variable  x that  measures  a 
characteristic  in  a finite  population  of  N elements.  Let 


(0.1) 


N p N p 

“ = if  J>i  * % = I I » 


be  the  population  mean  and  variance  of  the  measurements.  In  the  present 

2 

paper  the  problem  of  estimating  crN  on  the  basis  of  a sample  X^...^  , 

2 < n < N,  from  the  population  is  studied.  The  commonly  used  estimator 

is  the  sample  variance  crn  = ^ J.  (xi"xn)  > ^iere  xn  = ^ £ X..  is  the 

i=l  i=l  1 

sample  mean.  It  is  well  known  that  cr2  is  an  unbiased  estimator  of  cr2, 

under  simple  random  sampling. 

In  the  present  study  the  "unbiased"  estimator  <x  is  replaced  by 
a2  i n _ n 

^n  = n ^Xi“Xn^  * which  is  called  the  "classical"  estimator.  The  two 

estimators  are  nearly  equivalent  if  the  sample  size  is  not  trivially  small. 

The  "classical"  estimator  does  not  utilize  any  prior  information  on  or2  that 

may  be  often  available.  There  are  many  examples  of  repetitive  sampling 

surveys,  in  agricultural  or  industrial  areas,  in  which  good  information  is 

available  on  the  distribution  of  the  seasonal  or  yearly  yield  (production) 

of  a certain  commodity.  Samples  may  be  taken  during  the  season  to  observe 

the  distribution  of  related  characteristics  that  may  improve  the  forecasting 

of  a population  value.  Estimates  of  the  variance  in  the  population  could  be 

adjusted  adaptively  in  order  to  imporve  the  prediction  (confidence)  intervals 

for  population  parameters. 
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In  this  paper  we  show  how  such  prior  information  on  the  mean  and 
variance  of  the  population  can  be  utilized  to  adjust  the  "classical" 
estimator.  Specifically,  by  utilizing  the  special  structure  of  the  sample 
survey  theoretical  model  and  of  the  likelihood  function  we  derive  the 
general  form  of  Bayes  and  Bayes  Equi variant  estimators.  It  is  shown 
that  for  any  prior  distribution,  H,  of  independent  identically  distri- 
buted variables  x^, . . . ,x^,  having  a prior  mean  and  prior  variance 

cTq,  the  Bayes  estimator  of  cr  for  squared-error  loss  is  (approximately) 


(0.2) 


a2  n a2 
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Estimator  (0.2)  does  not  depend  on  the  sampling  procedure.  This  Bayes 
estimator  is  a wei^ited  average  of  the  "classical"  estimator,  o^,  based  on 
the  observed  sample  and  the  Bayes  estimator  of  the  "within"  variance  in  the 
unobserved  portion  of  the  population  and  the  estimator  of  the  variance 
between  the  means  of  the  observed  and  unobserved  portions  of  the  population. 

The  estimator  (0.2)  could  well  be  found  very  meaningful  and  good  also  in  a non- 

O 

Bayesian  sense  by  considering  and  Oq  as  proper  estimates  (or  quesses)  of 

the  mean  and  variance  of  the  unobserved  part  of  the  population. 

2 

Equivariant  estimators  of  the  variance  cr^.  are  considered  with  respect 

to  the  group  of  real  affine  transformations  on  the  parameter  space  9f^ 

2 

of  (x1, . . . ,Xjj).  It  is  shown  that  every  equivariant  estimator  of  cr^  can  be 

>\2 

expressed  in  the  general  form  (£n ) > where  \|r(un)  is  a proper  function 

of  the  maximal  invariant  statistic,  which  is  the  vector  of  standardized 

sample  values,  Bayes  equivariant  estimators  are  studied,  with  respect  to 

the  quadratic  loss  function  L(cr  \|r(u  )cr,T)  = (a  \|r(u  )-cr1.T)  /cr  . In  contrast 

v nT  v~n/  XT  n ~n'  XT  n 
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to  the  case  of  determining  Bayes  estimators,  the  form  of  the  Bayes  equivariant 
estimator  depends  strongly  on  the  particular  prior  distribution  specified 
for  x^,...,x^.  For  example,  it  is  shown  that 


(0.3) 


a2 
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n-3 
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a2  n 
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is  the  Bayes  equivariant  estimator  for  prior  normal  i.i.d.  variables, 
regardless  of  their  prior  mean  and  variance.  The  above  formula  (0.3 ) 
is  relatively  simple.  It  depends  only  on  ?rn  and  does  not  depend  on 
un.  This  is  not  always  the  case,  as  shown  by  Zacks  [16]  in  the  case  of 
exponentially  distributed  i.i.d.  variates. 

In  Section  1 we  introduce  the  sample  survey  model  and  discuss  sample 
statistics  and  likelihood  functions.  The  Bayesian  approach  extending  the 
sample  survey  model  is  discussed  in  Section  2.  Sections  3 and  4 define  and 
analyze  equivariant  estimators,  Bayesian  concepts  and  loss  functions. 
Bayesian  measures  of  relative  efficiency  are  introduced  in  Section  5. 

We  provide  a numerical  example  in  which  fifty  populations  of  size  N = 100 
were  simulated  from  an  exponential  distribution.  From  each  such  popula- 

a 2 2 

tion  a sample  of  size  n = 10  was  drawn  and  the  estimators  cr  and  aZ 

n B 

were  computed.  It  is  interesting  to  observe  the  extent  to  which  the 

Bayes  estimator  ou  is  more  effective  than  the  classical  estimator  cr 

B n 

in  small  samples.  Estimate  of  their  Bayes  relative  efficiency  is  provided 
in  that  example.  General  efficiency  analysis  is  provided  for  prior  normal 
distributions.  It  is  shown  that  the  classical  estimator  is  considerably 
less  efficient  than  the  Bayes  estimator.  Some  sensitivity  analysis  is 
performed  to  study  the  effects  of  erroneous  prior  parameters  on  the 
relative  efficiency. 
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There  are  only  a few  published  papers  on  the  Bayesian  estimation  of  the 

variance  of  a finite  population.  Liu  [9]  considered  unbiased  estimators 
2 

of  crw  under  various  possible  sampling  designs.  Since  the  sampling 
variance  of  these  estimators  depends  on  the  population  values  x1,...,x  , 

Liu  considered  the  Bayes  risk  of  these  estimators.  He  derived  a lower 
bound  to  the  Bayes  risk  function  and  thus  showed  the  optimality  of  the 
Horwitz -Thompson  type  estimator,  under  certain  conditions.  We  remark 
that  Liu's  study  is  not  really  a Bayesian  study,  since  proper  Bayes 
estimators  are  independent  of  the  sampling  design  and  are  generally 
not  unbiased.  Royall  [10],  [11]  and  Royall  and  Cumberland  [12]  studied 
the  problem  of  developing  confidence  intervals  of  the  population 
variance  by  regression  estimates.  We  cannot  compare  their  results 
with  ours  since  the  problems  are  different  and  so  are  the  approaches. 


2.  Foundations. 

Consider  a finite  population  of  size  U whose  units  have  values 

(real  finite)  x^,...,x^.  According  to  the  modern,  theory  of  sampling 

surveys  (see  Godambe  [5,6, 7]>  Basu  [2]  and  others)  the  population  vector 

)$H  = • • • >x-^)  considered  a parametric  point  in  a parameter  space 

(N) 

, which  belongs  to  the  Euclidean  N-space.  In  the  present  paper  a 
sample,  s,  of  size  n,  1 < n < N,  designates  a subvector  of  x„T 
consisting  of  n components  s = < x±  , . . . ,x.  >,  where 

€ [1^2,..., N}  for  all  j = 1, ...,N.  A sampling  procedure  is  a plan 
according  to  which  the  components  of  x^.  are  chosen.  In  a non-Bayesian 
theory  of  sampling  surveys  one  has  to  introduce  probability  functions 


4 


P(s)  on  the  sample  space,  %,  of  all  possible  samples,  in  order  to 
discuss  random  samples.  In  a Bayesian  theory  the  parametric  vector 
is  considered  a random  vector  having  a prior  joint  distribution 
According  to  this  approach,  the  population  vector. 


H(x,t)  on  % 


(N) 


Xjj,  is  a realization  of  a sample  from  a " superpopulation",  generated 

(like  in  a Monte  Carlo  procedure)  according  to  H(x^).  According  to 

this  approach,  given  any  sample  s = < x.  ,...,x.  >,  the  joint  prior 

1 1n 

distribution  of  s can  be  derived  from  H(x  ) and  the  posterior  joint 
distribution  of  xjj_n  = < Xy*  V ^ independent  of  the  sampling 

probability  function  P(s),  which  is  immaterial  for  a Bayesian  analysis 
(see  Solomon  and  Zacks  (1970)).  For  this  reason  we  will  assume  in  what 
follows,  without  loss  of  generality,  that  the  sample  consists  of  the 

>xj 

are  assumed  to  be  priorly  independent  and  identically  distributed  then 
any  sample  s can  be  considered  a simple  random  sample  from  H,  as 
in  the  classical  model  of  inference. 

The  estimation  problem  is  that  of  estimating  a specified  parametric 
function  e(x^)  of  the  population  vector  ( e.g.  the  population  mean, 
variance,  etc. ). 


sub  vector  Jn  = (x^ . . . ,xn)  and  ^.n  = <Vl’  • • • .*»>*  If 


3.  Estimators  of  the  Population  Variance. 

3*1  General  Structure. 

Let  £n  = • • • .>xn)  be  an  observed  sample.  Designate  by  x^, 

the  sample  mean  and  the  sample  variance,  respectively;  where 


- = 1 y 

n_ni=l  1 


and 


a2  1 v t - n2 

°n  - ; .1  W • 


i=l 
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If  x^_n  and  T^_n  designate  the  mean  and  variance  of  the  fen  units 
not  in  the  sample  then  the  population  variance  is 


(3.D  o^) . § • **  ♦ a- ♦ 1 a - 


Formula  (3«l)  can  be  verified  since  total  variance  may  be  written  as 
the  average  of  the  conditional  variances  plus  the  variance  of  the  condi- 
tional' expectations. 


Estimators  of  the  population  variance  are  sample  statistics  with  -range 

in  (0  ,oo).  The  most  common  estimators  in  use  are  the  sample  variance 

01  or  the  "unbiased  estimator"  o 2 - -—-a-2.  Formula  (3.1) 
n n n-1  n v ' 

shows- that,  regardless  of  the  sampling  procedure,  a proper  estimator  of 


^ (x^)  can  be  obtained  by  substituting  estimators  of  t2  anri 


— * 

fen  Xfen  in 

A a 

(3.1).  The  "unbiased"  estimator  can  be  obtained  from  (3.1)  by 

2 n 

fen  = HIT  • T^r  • 'an  • ^ wil1  be  shown 


substituting  xj_n  = xn  and  = ~ . ^-t1  . £2 


in  Section  4,  Bayes  estimators  of  a (x„)  can  be  obtained 

by  substituting  corresponding  Bayes  estimators  for  t2  and  x* 

fen  fen 

in  (3.1). 


3.2  Equi  variant  Estimators. 

Following  Fraser  [4]  we  will  denote  by  [a,p],  with  -co  < a < oo 

and  P ^ 0,  a real  affine  transformation,  i.e.  [a,P]x  = a+Px.  Let 

denote  the  group  of  all  such  transformations.  We  define  [a,P]x,T=(v-, ). 

where  y^=  [a,p]x^,  i=l, . . . , W.  Every  element  of  transforms  ST^  into 

in  a 1:1  fashion.  Let  & be  the  group  of  transformations  on  ' the  parameter 
2 

space  of  cr  (x^.)  induced  by  the  elements  of  Jf.  That  is,  if 
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kd  2 2 ? 

[a,p]x^  = yw  then  cr  (yN)  = P cr  (x^)  where  £ is  the  element  of^  corres- 

2 

ponding  to  [a,P]  of  G.  An  estimator  cr  (xn)  is  called  equi variant 
with  respect  to  G if,  for  every  [a,  [3]  e 


(3.2) 


a2([a,P]xn)  = P2o-2(xn)  , xn 


e 9T 


(n) 


The  sample  variance  an  is  equivariant  with  respect  to  )£.  The 
statistic 


x 


~n 


= [■ 


n 


cr 


n 


is  maximal  invariant  with  respect  to  jt^.  Thus,  every  equivariant  estimator 

2 

of  cr  (x^)  can  be  expressed  in  the  form 


(3.3) 


where  t(un)  is  a proper  positive  function  of  the  maximal  invariant 

statistic  u . For  further  reading  on  invariance  structures  for 
~n 

sampling  from  finite  populations  see  Chaudhuri  [3]. 


4.  Bayes  and  Bayes  Equivariant  Estimators. 

4.1  Bayes  Estimators. 

Let  H (x^)  be  a prior  distribution  in  a specified  family  ft  „ Let 

2 2 2 
L (cr  ,cr  ) denote  a loss  function  associated  with  estimating  cr  (x^)  by 

u2(x  ).  An  estimator  o-2(x  ) is  Bayes  with  respect  to  H(x.t)  and 
2 2 

L(cr  , cr‘  ) if  it  minimizes  the  prior  risk  function 


(4.1) 


R(cr  ,H)  = 


f 

J9f 


(N) 


L (a2  (x  ) 
v v~n 


2/ 

,cr  (x 


•N 


))dH(xK) 
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The  following  is  a general  result  for  the  squared-error  loss 


function: 


If  X.^,  ...  jXjj 
distribution  H(x) 
Bayes  estimator  of 


are  i.i.d.  random  variables  having  any  prior 

2 

with  a finite  prior  variance,  cr0,  then  the 
0-2  (Sfj)  is 


(4.2) 


A2  _ n A2  / n>r  2 n 
~ w + N (xn 


B N n 


dJ2] 


2 2 1 

viiere  cr^  = cTq  (l  - — ) and  is  the  prior  expectation. 

The  proof  of  (4.2)  proceeds  as  follows.  The  Bayes  estimator  of 

2 

17  §iven  for  the  squared-error  loss  function  is  the  posterior 

expectation  of  (3.1).  Furthermore,  since  the  components  of  x,  are 
priorly  independent 


(4.3) 


~n^  'VTN-n^  " 


H-n-1  2 

N-n  ^0 
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for  any  prior  distribution  H,  having  variance  crn.  Moreover, 


(4.4) 


Substituting  these  expressions  in  (3.1)  one  obtains  (4.2).  In 
many  situations  it  is  not  unreasonable  to  assume  that  x^, . . . ,x^  are 
priorly  i.i.d.  Hence,  formula  (4.2)  is  a very  general  formula,  since 
it  does  not  depend  on  the  form  of  H(x),  but  only  on  the  prior  mean 
and  variance.  These  values  may  be  known  from  previous  experience. 
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4.2  Bayes  Equi variant  Estimators. 

Consider  the  structure  of  Bayes  equi variant  estimators.  We  have 

* r Xn  I-,  * 

-Sn-ii  a > a J£N_n  > 

n n 

where  x^_n  = (xn+1>  • • • jX^)  and  is  maximal  invariant  with  respect  to 
Mr>  Let  vN_n  and  w2_n  be  the  mean  and  the  variance  of  u*  One  can 
express  the  population  variance  indieae  terms  in  the  form 


(4.  5) 


2;  % A2rn  , 

(2n>  ’ °n[it  + & - 


nw  2 
^WN-n 


, n 2 ,, 

N N-n'J 


Thus,  comparing  (3.3)  and  (4.5),  the  t-fnnction  of  an  equi variant 
estimator  should  be  chosen  to  estimate  the  function 


(4.6) 


D(w 


Nrn 


2 , _ n 

,vN-n;  ~ I 


+ 0<"L 


, n 2 \ 

IT  vN-n^ 


Let  L(t,D)  be  a loss  function  for  the  estimation  of  D(w2  v2  ) bv 

N-n*  N~ny  y 

+ • . L(feD)  is  invariant  with  respect  to  Jf  . Let  G(un,u|  ) be 

a prior  distribution  induced  by  H(^).  ae  prior  risk  associated  with 
t and  G is 


(4.7) 


R(t,G)  = JWi(un);D(w|_nJv|_n))aG(Vuj_n; 


2 

An  estimator  crntG(un)  is  called  Bayes  equivariant  if  minimizes  (4.7). 

Notice  that  the  criterion  of  minimizing  (4.7)  is  the  same  as  minimizing 

the  Bayes  risk  for  the  quadratic  loss  L(e,cr2T)  = (£-cr2)  /a,  where 

N N n 

® = V^)* 
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In  many  applications  it  would  be  reasonable  to  assume  that  the  family 
of  prior  distributions  is  a family  with  location  and  scale  parameters.  In 
other  words,  assume  that  all  the  prior  distributions  of  are  of  the  form 


H( 


X1^0  X2^0 


0 


0 


) , 


where  -oo  < ^ < oo  and  0 < cr^  < co.  In  this  case  the  Bayes  equi variant 

estimator  depends  only  on  the  general  form  of  H(x^, . . . ,x^) . Indeed*.-;  the 
distribution  G(un,u^  n)  is  the  same  for  all  and  cr^  of  distributions 

in  Hr. 

Equivariant  estimators  in  the  strict  sense  were  defined  as  those  of 
the  form  or  \|r(u  ).  The  Bayes  estimator  cr^  (4.2)  is  thus  not  strictly 
equivariant.  However,  if  x is  transformed  to  [o^PJx^  the  prior 

parameters  (m-q-*0^  should  be  transformed  to  [«,P]  (p0>°-q ) = (a+Pq,-,*  I P I crQ). 

A 2 A 2 

Let  cr  / s (x  ) denote  the  Bayes  estimator  cr^  with  the  prior  parameters 

(mo^O  ~n  B 

w w a2 

q^  and  crn,  respectively.  Then  a , _ N (x_ ) is  generalized  equivariant  in 

u u ~n 

the  sense  that 


(4.8) 


a?  RU  n([«,P]x  ) = P2&^  Ax  ) , 
[a*PJ  (.M-q^o)  (,4n*crn;  ~n 


0' 


for  all  *oo  < a < oo,  0 < P < oo;  and  all  xn<  Furthermore,  the  Bayes 


estimator  ov  \ v.x. 

(q  ,a  )v~n 

u u * 2 \2  /A  4 

equivariant  estimators  with  respect  to  the  quadratic  loss  (e-cxTO)  /°"n' 


(x^)  is  also  Bayes  in  the  class  of  all  generalized 

2 > 
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4.3  Examples  of  Bayes  Equi variant  Estimators. 


4.3.1  Normal  Priors. 


Suppose  that  are  priorly  independent  and  identically 

distributed  (i.i.d. ) normal  variables  with  prior  mean  p and  prior 

2 ^ 
variance  cr0.  % consists  of  all  such  distributions  with  -00  < p < 
.2  0 
and  0 < o-q  < 00  . Under  this  model  the  sample  statistics  x and  cr2 

n n 


CO 


are  independent  of  t 

N»n 


and  . The  Bayes  equivariant  estimation  is 


(4.9) 


°BE  ~ °"n{N  + (l  - |) [E{i w§_Jun}  + I E{vN-J)2n)]  3 


We  now  show  that  w|_n  and  v^_n  are  independent  of  Indeed,  by 

the  Bayes  model  xn  and  x^_n  are  independent.  Hence  &n>vn,\)  is 

independent  of  x^_n-  Furthermore,  (xn,£n)  is  a complete  sufficient 

statistic  for  the  subfamily  of  prior  distributions  of  x . Hence,  from 

~n 

Basu’s  theorem  (Basu,  [l] ) un  is  independent  of  (x^,!^).  Finally, 

since  Bw-n  is  a function  of  fe^n'V^n^  £n  and  inde- 

pendent. Hence, 


(4.10) 


E{F[N-n-l,n-l]J 


n . N-n-1 
n-3  N-n  * 


Similarly, 
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(4.ii) 


E^VN-J~n}  EfvN-n5  (N-n)(n-3) 


Substituting  these  results  in  (4.9)  we  obtain  as  the  Bayes  equi variant 
estimator 


(4.12) 


tW$n^  " N^1  + n-3  ^ 9 


2 

which,  in  large  populations  is  close  to  En  (x--x)  /(n-3). 

i=l  1 

It  is  well  known  (see  Zacks  [17 J pp.  346))  that  the  mi ni m-nm  mean- 

2 

squared-error  equi variant  estimator  of  cr^  in  the  i.i.d.  case  is 
n wmmm  2 

1 . (x^-x)  /(n+l).  On  the  other  hand,  if  the  loss  function  is  the 

A 2 2 A U 

quadratic  loss  (6-cTq)  /cr^  the  best  equi  variant  estimator  is 

■Q  mb  2 

li=1  (x^x)  /(n-3).  This  confirms  the  above  result. 

4.3.2  Exponential  Priors 

Suppose  that  x^, . . . ,x^  are  priorly  i.i.d.,  with  a common  exponen- 

2 2 

tial  distribution,  with  mean  ^ (the  prior  variance  is  cr0  = pQ).  It 
is  shown  in  [l6]  that  the  Bayes  equi variant  estimator  is 


(4.13) 


f2  (u  )n2 
(n-2 ) (n-3  ) 


> 


where  f(un)  = -min(u^,  ...,u  ).  Thus,  in  contrast  to  the  normal  case, 

in  the  exponential  priors  model  the  Bayes  equi variant  estimator  depends 

on  u too. 

~n 
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5*  Bayes  Risk  Efficiency. 

In  the  present  section  we  introduce  an  index  of  Bayesian  efficiency 

of  estimators  of  cr^.  Given  a prior  distribution,  HCx^.),  we  denote  by 

R(cr  ,H)  the  -prior  risk  function  of  an  estimator  a2  under  H.  In  the 

present  section  we  consider  squared-error  loss,  (o-2-cr2)2,  only. 

Generalization  to  quadratic  loss  functions  can  be  readily  attained. 

Since  the  minimal  prior  risk  is  attained  by  the  Bayes  estimator  a2, 

B 

o 

with  proper  prior  parameters  p , cr^,  we  define  the  prior 
relative  efficiency  of  an  estimator  cr2  as 


(5.1) 


ke(£2,h) 


R(^H) 

r(£2,h) 


For  any  estimator,  0 < RE(cr2,H)  < 1. 


5*1  Exponential  Priors -Example. 

We  provide  now  a numerical  example  of  estimating  the  variance  of 

a small  population,  N = 100,  When  the  population  variates  are  priorly 

i.i.d.  exponential  random  variables  with  expectation  p^  = 10.  The  sample 

size  is  n = 10.  In  Table  1 we  present  the  values  of  the  classical  and 

the  Bayes  estimators  determined  by  50  independent  simulation  runs.  In 

each  case  we  give  also  the  value  of  cr|.  We  see  that  generally  the 

Bayes  estimator  is  closer  to  the  population  variance.  The  prior  relative 

efficiency  of  o'2  against  o'2  is  estimated  to  be  RE  = .122.  Note  R(cr2) 
n B n 

a /*2> 

and  RpjgJ  are  the  sample  estimates  of  the  prior  mean-squared-errors 
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a2  2 2 a2  2 2 

V (o-n-^N)  ) and  EfjUo-g-OTjj)  h respectively.  We  see  in  this  example 
that  the  classical  sample  variance  is  very  inefficient  compared  to  the 
Bayes  estimator.  In  the  following  example  we  show  some  analytical 
comparisons  for  the  normal  case. 


5*2  Normal  Priors . 

The  prior  relative  efficiency  index  (5*l)  can  be  expressed  also 
in  the  form 


(5-2) 


A2  y (£2-£|)2}  _i 

>E)  = [1  + \lPVBj  3 


where  PVR  is  the  posterior  variance  of  the  Bayes  estimator  cr2. 

B 

the  case  of  prior  i.i.d.  normal  (M-n>  cr2)  variables , one  obtains 
(5-3) 


In 


2a. 


'O'  Cr 

4 
0 


E{PVE(H,xn)}  = (l-f)(l-i(l-f))  , 


where  f = n/u.  Consider  the  sample  variance  a2.  Due  to  the  prior 
“■  2 

independence  of  x and  cr  we  obtain 

n n 

(w)  W” 


2 2 

where  X^[n-1]  and  X^[l]  designate  independent  chi -squared  r.v. 's. 

From  (5.3)  and  (5.4),  the  prior  relative  efficiency  of  cr2,  relative 

n 

14 


Table  1.  50  Independent  Simulation  Runs  of  Exponential  Populations 

of  Size  N = 100.  Variance  Estimates  are  Based  on  Sample 
of  Size  n = 10. 


A2f-\ 

cr  (1 ) 


A2  /.  s 
^00 


2,.  \ 


1 

11.0279 

91.3822 

2 

25.9047 

93.5383 

3 

163.1974 

108.4712 

4 

114.887 6 

102.4455 

5 

48. 5473 

96.1006 

6 

24.1673 

92.6005 

7 

233.9413 

115.0390 

8 

24.2149 

92.8096 

9 

36.0928 

93.7894 

10 

58.1589 

95.8768 

ll 

82.4155 

98.2541 

12 

110.1057 

101.4133 

13 

30.7288 

94.2891 

14 

92.7353 

99.3392 

15 

14.1656 

92.9971 

16 

30.6368 

93.1324 

17 

75.4294 

97.7225 

18 

123.1242 

102.-9664 

19 

34.9740 

95.1125 

20 

226.5559 

116.6118 

21 

16.9147 

92.2562 

22 

98.5459 

100.9711 

23 

48.8585 

94. 9472 

24 

54.3939 

95.4909 

25 

52.1598 

95-2457 

26 

55.5302 

95.6501 

27 

171.1030 

107.1958 

28 

105.0818 

100.6807 

29 

17.5409 

92.9367 

30 

23.1610 

93.2519 

31 

116.6341 

101.8029 

32 

27.0976 

94.1436 

33 

69.4181 

98.1458 

34 

20.7325 

93.0247 

35 

191.2939 

112. 5688 

36 

199. 9331 

110.6755 

37 

83.3183 

98.5574 

38 

17.2250 

95.0543 

39 

218.7134 

112.1939 

4o 

30.0034 

93-1121 

4l 

34.1185 

93.4987 

42 

32.7168 

93.6423 

43 

50.9459 

95.3254 

44 

41.4210 

94.8028 

45 

36.6681 

94.1386 

46 

22.9323 

93 . 6664 

47 

48 

41.4523 

149.5766 

94.1452 

106.0161 

49 

19.2172 

101.9416 

50 

64.9260 

96.6517 

59-3401 

99-7939 

117-9033 

112.7694 

140.4537 
95-3784 

110.2580 

73.1228 

88.0842 

92.6123 

100.8301 

82.7390 

82.7665 

66.7146 

100.5613 

87.0248 
74.1037 
104. 8201 

86.4942 

118.4538 
134.9088 

72.2479 

134.9985 

ll4.8o43 

158.8236 

79.0303 

87.0796 
128.0554 
95.6511 
94.2407 
118.2461 
74.1948 
101.3439 
98.1385 
121.3687 
105.3757 
67 . 86i4 
77.5798 
108.5281 
83.3298 
80.7266 
106.1845 
77.0860 
64.2112 
75.0202 
88. 6936 
64.5216 
99.7121 
79.8605 
60.5816 


R(a^)  = 3887.67 

R(c|)  = 474.90 
RE(a^H)  = .12216 
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to  the  prior  normal  distributions  is: 


(5-5) 


ke(&2,h) 


1 + 


N(l-f)[-|  + (1-|)2  - 
2(1  - |(l-f)) 


The  relative  efficiency  function  is  independent  of  the  prior 
parameters,  since  o-Q  is  a scale  parameter  of  the  distribution.  We 
therefore  provide  in  the  following  table  some  relative  efficiency  values 
as  functions  of  the  sample  fraction  f and  the  population  size,  N. 


Table  2 . The  Prior  Relative  Efficiency  of  cr  . 


N\f 

0.10 

0.25 

0.50 

0.75 

100. 

o.n4 

0.310 

0.666 

0.921 

200. 

0.112 

0.309 

0.666 

0.922 

300. 

0.111 

0.308 

0.666 

0.922 

4oo. 

0.111 

0.308 

0.666 

0.922 

500. 

0.111 

0.308 

0.666 

0.923 

600. 

0.111 

0.308 

0. 666 

0.923 

700. 

0.111 

0.308 

0.667 

0.923 

800. 

0.110 

0.308 

0.667 

0.923 

900. 

0.110 

0.308 

0.667 

0.923 

1000. 

0.110 

0.308 

0.667 

0.923 

From  Table  2 the  prior  relative  efficiency  of  cr^  is  almost  independent  of 
the  population  size  N and  is  somewhat  greater  than  the  sample  fraction, 
f.  These  numerical  results  show  the  extent  of  possible  improvement  in 
estimation  if  good  information  is  available  on  the  prior  distribution. 


16 


In  order  to  analyze  the  extent  of  errors  in  the  prior  assumptions  concerning 

the  values  of  and  cr0  we  derive,  on  the  basis  of  (5.2),  the  prior 

a2 

relative  efficiencies  of  cr  (x  ),  under  (p~,ou).  It  is  a straight- 

p.,cr.  ~n  O 0 

^ A 2 

forward  matter  to  show  that  the  prior  relative  efficiency  of  cr  (x  ) 

lV°l  n 

is 

2 2 

P W(l-f)[(p-l+fB2)  -4f2  ~]  -i 

(5.6)  KE  (S’2  H)  = [1  + = — ] , 

H'  2(1-|  (1-f)) 

2 2 

where  p = cr^/c^  and  6 = (p1-pQ)/cr0. 

In  Table  3 we  present  the  prior  relative  efficiency  of  the  Bayes  estimator 

2 2 
cr^,  (xn)  as  a function  of  f,  & and  X=p-1,  where  H*  is  the  N(jjt^,cr^)  distribution. 

We  see  that  the  magnitude  ©f  & is  not  so  important,  but  deviations  from  cr0 

larger  in  magnitude  than  10  percent  reduce  the  prior  relative  efficiency 

a 2 

below  that  of  crn.  In  Table  4 we  provide  these  prior  relative  efficiency 
values  for  values  of  \ between  -7-5$  to  7.5$  • We  see  that  in  this 
range  the  Bayes  estimator  is  considerably  more  efficient  than  the  classical 
sample  variance. 
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Table  3.  The  Frior  Relative  Efficiency  of  *h«  (xn)>  N=  1,000 


Sample  Fraction  = 

= .10 

b\x 

-.30 

— # 20 

-.10 

0.00 

0.10 

0.20 

0.30 

-.50 

0.028 

0.067 

0.279 

0.754 

0.124 

0.042 

0.021 

-.1*0 

0.027 

0.061 

0.238 

0.874 

o.i4i 

0.045 

0.022 

-.30 

0.026 

0.057 

0.211 

0.950 

0.157 

0.048 

0.023 

-.20 

0.025 

0.055 

0.194 

0.986 

0.170 

0.051 

0.023 

-.10 

0.021+ 

0.053 

O.I85 

0.998 

0.179 

0.052 

0.024 

.00 

0.021+ 

0.053 

0.182 

1.000 

0.182 

0.053 

0.024 

0.10 

0.021+ 

0.053 

0.135 

0.998 

0.179 

0.052 

0.024 

0.20 

0.025 

0.055 

0.194 

0.986 

0.170 

0.051 

0.023 

0.30 

0.026 

0.057 

0.211 

0.950 

0.157 

0.048 

0.023 

0.1+0 

0.027 

0.061 

0.238 

0.874 

o.i4i 

0.045 

0.022 

0.50 

0.028 

0.067 

0.279 

0.754 

0.124 

0.042 

0.021 

Sample  Fraction  = 

= .25 

&\X 

-.30 

-.20 

-.10 

0.00 

0.10 

0.20 

0.30 

-.50 

0.01+5 

0.122 

0.617 

0.391 

0.091 

0.037 

0.020 

-.1+0 

0.038 

0.094 

0.415 

0.602 

0.119 

0.044 

0.023 

-.30 

0.033 

0.078 

0.304 

0.817 

0.150 

0.051 

0.025 

-.20 

0.031 

0.069 

0.247 

O.950 

0.180 

0.057 

0.027 

-.10 

0.029 

o.o64 

0.219 

0.99*4- 

0.202 

0.061 

0.028 

.00 

0.029 

0.062 

0.210 

1.000 

0.210 

0.062 

0.029 

0.10 

0.029 

o.o64 

0.219 

0.994 

0.202 

0.061 

0.028 

0.20 

0.031 

0.069 

0.247 

0.950 

0.180 

0.057 

0.027 

0.30 

0.033 

0.078 

0.304 

0.817 

0.150 

0.051 

0.025 

0.1+0 

0.038 

0.094 

0.415 

0.602 

0.119 

0.044 

0.023 

0.50 

0.01+5 

0.122 

0.617 

0.391 

0.091 

0.037 

0.020 

Sample  Fraction  = 

= .50 

5\X 

-.30 

-.20 

-.10 

0.00 

0.10 

0.20 

0.30 

-.50 

o.ii4 

0.395 

0.780 

0.199 

0.073 

0.036 

0.022 

-.1+0 

0.076 

0.214 

0.847 

0.373 

0.109 

0.048 

0.027 

-.30 

0.058 

0.142 

0.555 

0.645 

0.159 

0.062 

0.032 

-.20 

o.o48 

0.110 

0.382 

0.893 

0.216 

0.076 

0.038 

-.10 

o.o44 

0.095 

0.307 

0.989 

0.266 

0.087 

o.o4i 

.00 

0.043 

0.091 

0.286 

1.000 

0.286 

0.091 

0.043 

0.10 

0.044 

0.095 

0.307 

0.989 

0.266 

0.087 

o.o4i 

0.20 

o.o48 

0.110 

0.382 

0.893 

0.216 

0.076 

0.038 

0.30 

0.058 

0.142 

0.555 

0.645 

0.159 

0.062 

0.032 

0.1+0 

0.076 

0.214 

0.847 

0.373 

0.109 

o.o48 

0.027 

0.50 

0.114 

0.395 

0.780 

0.199 

0.073 

0.036 

0.022 
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Table  4.  Prior  Relative  Efficiency  of  ct^,  (x^)  , E=  1,000. 


Sample  Fraction  = 

.10 

&\x 

-.075 

-.050 

-.025 

0.000 

0.025 

0.050 

0.075 

-.50 

o.46i 

0.754 

0.957 

0.754 

o.46i 

0.279 

0.180 

-.40 

0.385 

0.645 

0.939 

0.874 

0.560 

0.334 

0.210 

-.30 

0.336 

0.564 

0.884 

0.950 

0.651 

0.387 

0.238 

-.20 

0.305 

0.510 

0.829 

0.986 

0.722 

0.431 

0.262 

-.10 

0.288 

0.480 

0.793 

0.998 

0.766 

o.46o 

0.278 

-.00 

0.283 

0.470 

0.780 

1.000 

0.780 

0.470 

0.283 

0.10 

0.288 

0.480 

0.793 

0.998 

0.766 

o.46o 

0.278 

0.20 

0.305 

0.510 

0.829 

0.986 

0.722 

0.431 

0.262 

0.30 

0.336 

0.564 

0.884 

0.950 

0.651 

0.387 

0.238 

o.4o 

0.385 

0.645 

0.939 

0.874 

0.560 

0.334 

0.210 

0.50 

0.46i 

0.754 

0.957 

0.754 

o.46i 

0.279 

0.180 

Sample  Fraction  = 

.25 

sU 

-.075 

-.050 

-.025 

0.000 

0.025 

0.050 

0.075 

-.50 

0.868 

0.868 

O.617 

0.391 

0.252 

0.171 

0.122 

-.1+0 

0.658 

0.911 

0.874 

0.602 

0.378 

0.244 

0.166 

-.30 

0.484 

0.759 

O.965 

0.817 

0.532 

0.333 

0.217 

-.20 

0.385 

0.619 

0.910 

0.950 

0.678 

0.423 

0.268 

-.10 

0.336 

0.540 

0.838 

0.994 

0.777 

0.491 

0.307 

-.00 

0.321 

0.516 

0.810 

1.000 

0.810 

0.516 

0.321 

0.10 

0.336 

0.540 

0. 888 

0.994 

0.777 

0.491 

0.307 

0.20 

0.385 

0.619 

0.910 

0.950 

0.678 

0.423 

0.268 

0.30 

0.484 

0.759 

0.965 

0.817 

0.523 

0.333 

0.217 

o.4o 

0.658 

0.911 

0.874 

0.602 

0.378 

0.244 

0.166 

0.50 

0.868 

0.868 

0.617 

0.391 

0.252 

0.171 

0.122 

Sample  Fraction  = 

.50 

&u 

-.075 

-.050 

0.025 

0.000 

0.025 

0.050 

0.075 

-.50 

0.571 

0.395 

0.276 

0.199 

0.148 

0.114 

0.090 

-.40 

0.921 

0.766 

0.544 

0.373 

0.261 

0.188 

o.i4l 

-.30 

0.787 

0.951 

0.873 

0.645 

0.440 

0.303 

0.215 

-.20 

0.573 

0.803 

0.974 

0.893 

0.655 

0.445 

0.305 

-.10 

0.448 

0.662 

O.905 

0.989 

0.813 

0.568 

0.384 

-*KX) 

.0.415 

O.615 

0.865 

1.000 

0.863 

0.615 

0.415 

0.10 

0.448 

0.662 

0.905 

0.989 

0.813 

0.568 

0.384 

0.20 

0.563 

0.803 

0.974 

0.893 

0.655 

0.445 

0.305 

0.30 

0.787 

0.951 

0.873 

0.645 

o.44o 

0.303 

0.214 

o.4o 

0.921 

O.766 

0.544 

0.373 

0.261 

0.188 

0.l4l 

0.50 

0.571 

0.395 

0.276 

0.199 

o.i48 

0.114 

0.090 
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